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ABSTRACT 


This  report  summarizes  a  series  of  tests  designed  to  determine  the  response  of 
quasi-dc  currents  on  three-phase  power  distribution  transformers  for  electric  power 
systems.  In  general,  if  the  dc  injection  is  limited  to  the  primary  side  of  a  step-down 
transformer,  significant  harmonic  distortion  is  noted  and  an  increase  in  the  reactive  power 
demand  results.  For  dc  injection  on  the  secondary  (load)  side  of  the  step-down  transformer, 
the  harmonic  content  at  the  secondary  side  is  quite  high  and  saturation  occurs  with  a 
relatively  low  level  of  dc  injection;  however,  the  reactive  power  demand  is  significantly 
lower.  These  tests  produced  no  apparent  damage  to  the  transformers.  Transformer  damage 
is  dependent  on  the  duration  of  the  dc  excitation,  the  level  of  the  excitation,  and  on  the 
thermal  characteristics  of  the  transformer. 

The  transformer  response  time  is  found  to  be  much  shorter  than  seen  in  power 
transformer  tests  at  lower  dc  injection  levels.  This  shorter  response  time  suggests  that  the 
response  time  is  strongly  dependent  on  the  injected  current  levels,  and  that  higher  levels  of 
dc  injection  for  shorter  durations  could  produce  very  high  reactive  power  demands  and 
harmonic  distortion  within  a  few  tenths  of  a  second. 

The  added  reactive  power  load  could  result  in  the  blowing  of  fuses  on  the  primary 
side  of  the  transformer  for  even  moderate  dc  injection  levels,  and  neutral  currents  are  quite 
large  under  even  low-level  dc  injection.  This  "smoking  neutral"  results  in  high-level 
harmonic  injection  into  equipment  via  the  neutral  and  in  possible  equipment  failure. 
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IMPACT  OF  QUASI-DC  CURRENTS 
ON  THREE-PHASE  DISTRIBUTION  TRANSFORMER  INSTALLATIONS 


1.0  OBJECTIVES 

It  is  a  well-established  fact  that  geomagnetic  storms  influence  electrical  power 
transmission  and  distribution  systems  [1-3].  Previous  cases  of  such  storms  in  the  northern 
latitudes  have  resulted  in  occasional  power  disruptions  and,  in  some  cases,  damage  to 
transformers  [4],  These  effects  are  caused  by  a  time  variation  of  the  earth's  magnetic  field 
creating  an  induced  electric  field  (E-field)  along  the  surface  of  the  earth.  This  E-field  acts  as 
a  voltage  source  along  long  power  transmission  or  distribution  lines,  and  if  such  lines  are 
connected  to  the  earth  at  both  ends,  a  quasi-static  current  can  flow.  This  current  can  cause 
unwanted  saturation  in  the  magnetic  cores  of  transformers  in  the  power  system,  which  in  turn 
produces  harmonic  distortion  and  transformer  heating.  This  harmonic  distortion  and  heating 
can  lead  to  system  upset  (shutdown)  and  possibly  transformer  burn-out. 

A  high-altitude  nuclear  explosion  is  also  known  to  affect  the  magnetosphere, 
producing  late-time  variations  of  the  earth's  magnetic  field  for  several  hundreds  of  seconds. 
Known  as  the  magnetohydrodynamic  electromagnetic  pulse  (MHD-EMP),  or  E3,  this 
magnetic  field  disruption  is  of  particular  concern  to  electrical  power  systems.  Its  effects  are 
similar  to  those  of  geomagnetic  storms.  Although  the  MHD-EMP-induced  currents  can  be 
significantly  larger  in  magnitude,  they  last  for  a  shorter  period  of  time  than  do  those  from  a 
geomagnetic  storm.  The  effects  of  the  MHD-EMP  environment  compound  the  adverse 
effects  of  the  early-time  high  altitude  EMP  (HEMP)  environment  [5],  posing  a  potentially 
serious  threat  to  the  electrical  power  system. 

In  order  to  increase  the  understanding  of  the  effects  of  the  MHD-EMP  environment 
on  electrical  power  systems,  several  investigations  have  been  conducted.  An  initial  study  [6] 
has  developed  the  framework  for  understanding  the  MHD-EMP  interaction  with  power 
systems,  and  outlines  a  methodology  for  assessing  the  system  behavior  in  this  environment. 
A  subsequent  report  [7]  continues  with  studies  of  MHD-EMP  coupling  to  realistic 
transmission  and  distribution  lines,  using  more  up-to-date  environments. 

The  present  report  documents  an  experimental  program  designed  to  increase  the 
understanding  of  the  behavior  of  distribution-class  power  transformers  subjected  to  quasi-dc 
current  excitation.  Given  the  knowledge  of  the  MHD-EMP-induced  current  flowing  in  a  long 
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power  line,  and  the  given  transformer  response  characteristics  determined  in  this  program, 
more  accurate  assessments  of  the  behavior  of  the  overall  power  system  to  EMP  can  be  made. 

The  experimental  program  was  conducted  in  June  of  1990  at  the  Mission  Research 
Corporation  in  Albuquerque,  NM,  and  involved  constructing  a  mock-up  of  a  simple  3-phase 
12.47-kV  distribution  system  using  two  distribution-class  power  transformers  and  a  dummy 
3-phase  load  bank.  This  system  was  fed  by  local  commercial  power.  Measurements  on  this 
system  were  made  with  the  following  specific  objectives  in  mind; 

1.  To  determine  the  effect  of  quasi-dc  currents  on  the  operation  of  three-phase  transformer 
banks. 

2.  To  measure  voltage  and  current  harmonics  within  the  system  and  at  the  loads. 

3.  To  assess  the  importance  of  the  quasi-dc  current  duration. 

4.  To  determine  the  change  in  reactive  power  demand  as  a  function  of  quasi-dc  currents. 

5.  To  determine  if  low-level  quasi-dc  currents  and  the  distorted  ac  current  can  cause  primary 
fuses  to  blow. 

Objectives  1,  2,  and  4  were  met  taking  data  directly  from  these  measurements,  while 
objectives  3  and  5  were  met  by  interpreting  the  data. 
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2.0  TEST  ORGANIZATION 


For  this  experimental  program,  the  Oak  Ridge  National  Laboratory  (ORNL)  was  the 
lead  organization.  ORNL  responsibilities  included  performing  the  functions  of  project 
coordination,  test  plan  development,  pretest  and  posttest  analysis,  and  writing  the  final  report. 
Assisting  ORNL  was  Dr.  F.  M.  Tesche,  a  private  consultant.  During  the  measurement  phase 
of  the  program,  support  was  provided  by  Aesa  Brown  Boveri  (ABB)  Power  Systems,  Inc. 
Mission  Research  Corporation  (MRC)  in  Albuquerque  provided  the  test  site,  and  all 
equipment  and  personnel  associated  with  the  data  measurement,  storage,  viewing,  and 
conversion.  MRC  also  provided  the  dc  current  injection  source  and  the  480-V  3-phase 
commercial  power  hook-up.  The  Public  Service  Company  of  New  Mexico  provided  the 
208-V/12.47-kV,  75-kVA  Y-Y  and  the  480-V/I2.47-kV,  300-kVA,  Y-Y  3-phase 
transformers  for  the  test.  In  addition,  MRC  obtained  a  12.47-kV/208V,  A-Y,  75-kVA  3- 
phase  transformer  under  a  purchase  order  fi-om  ABB. 


3.0  TEST  CONFIGURATION 

For  this  test,  a  mock-up  of  a  power  distribution  system  was  made  in  the  MRC  test  bay. 
An  existing  480-V  3-phase  commercial  power  service  was  available  for  the  primary  energy 
source.  This  energy  source  had  the  configuration  shown  in  Figure  1,  where  90-A  slow-blow 
fuses  served  as  primary  protection.  The  fuse  box  and  service  disconnect  switch  were 
connected  to  a  30  ft  1/0  cable  leading  to  the  test  area,  where  a  200- A  3-phase  contactor  was 
used  to  energize  and  deenergize  the  transformers  being  tested. 

Three  different  transformer  configurations  were  used  in  this  test.  These  are  shown  in 
Figures  2,  3,  and  4.  The  first  configuration,  denoted  as  lA,  involved  the  480-V/12.47-kV, 
300-kVA,  Y-Y  transformer,  shown  as  Tl,  and  the  12.47-kV/208-V,  75-kVA,  Y-Y 
transformer,  shown  as  T2.  The  neutrals  of  the  secondary  (12.47  kV  side)  of  Tl  and  the 
primary  (12.47  kV  side)  of  T2  were  connected  to  the  dc  injection  source,  which  simulated  the 
MHD-EMP  excitation  of  a  long  distribution  line  terminated  by  grounded  Y-Y  transformers  at 
each  end.  Three  resistive  dummy  loads  comprised  a  nominal  60-kW  load  for  transformer  T2. 
This  load  was  estimated  to  be  about  80%  of  the  design  load  for  this  transformer.  The  total 
60-Hz  current  in  each  phase  as  supplied  by  the  480-V  source  was  estimated  to  be  80  to  88  A. 
Table  1  summarizes  the  transformer  and  load  bank  characteristics. 
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Existing  Power 
Service 


Three-Phase 
Equipment  in 
Test  Area 


Figure  1 .  Commercial  480-V  Source  Configuration. 


Figure  2.  Test  Configuration  1  A. 


12  VOC 
CONTROL 


DC  injection  circuit 


Figure  3.  Test  Configuration  IB, 


Figure  4.  Test  Configuration  2 


Table  1. 

Transformer  and  Load  Bank  Characteristics 


Equipment 

Model  No. 

Description 

T1  GE  Transformer 

K860313T72AA 

300-kVA  12.47-kV/480-V  Y-Y  3<t> 

T2  GE  Transformer 

M318597TFPA 

75-kVA  12.47-kV/208-V  Y-Y  30 

T3  Weaver  Transformer 

23490 

75-kVA  12.47-kV/208-V  A-Y  30 

Load  Banks 

Three  55-gallon  drums  with  water 
heater  elements  in  parallel.  Dc 
resistance;  Phase  1  -  0.68  0,  Phase  2  | 
-  0.70  Q,  Phase  3  -  0.72  0. 

The  second  system  configuration  (IB)  shown  in  Figure  3,  involved  the  same 
transformers  as  configuration  lA,  but  the  quasi-dc  excitation  was  provided  in  the  secondary 
(208  V  side)  load  circuit  of  transformer  T2.  The  last  test  configuration,  denoted  by 
configuration  2  in  Figure  4,  shows  the  12.47-kV/208-V,  A-Y,  75-kVA  transformer  as  T3, 
with  the  quasi-dc  again  being  induced  in  the  load  side  of  T3.  In  each  of  these  configurations, 
the  primary  side  of  transformer  T1  is  considered  to  be  the  480  V  side,  and  the  secondary-side 
is  the  12.47  kV  side.  For  transformers  T2  and  T3,  the  primary-side  is  the  12.47  kV  section, 
and  the  secondary  is  the  208  V  side. 

For  the  purpose  of  possible  future  modeling  efforts,  we  wanted  to  obtain  a  suitable 
electrical  model  for  the  transformers  used  in  this  test.  Figure  Sa  illustrates  one  possible 
circuit  model  for  a  single  phase  of  the  three  phases  of  a  general  transformer  [8].  The  shunt 
resistance  and  inductance  elements  and  L,.  represent  the  core  loss  and  the  core  magnetizing 
inductance,  and  the  series  elements  and  represent  the  winding  resistance  and  the  leakage 
inductance  of  the  transformer  windings.  Connected  to  these  elements  is  an  ideal  transformer 
having  a  turns  ratio  of  l:a.  In  these  models,  we  assume  that  the  primary  is  the  low-voltage 
side  of  the  transformer  and  the  secondary  is  the  high-voltage  winding.  Thus,  the  turns  ratio 
"a"  is  greater  than  unity. 

At  the  conclusion  of  the  testing,  measurement  were  made  for  determining  the 
appropriate  parameters  for  the  three  transformers  used  in  the  experiment.  These  tests 
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involved  a  short-circuit  test  to  determine  and  Lj ,  as  well  as  an  open-circuit  test  for  R^.  and 
Lg.  These  tests  and  their  rationale  are  discussed  in  detail  in  ref.[8].  The  general  test  setup  for 
these  measurements  is  shown  in  Figures  5b  and  5c,  along  with  the  relations  for  determining 
the  parameters  G,  B,  R,  and  X.  These  parameters  are  related  to  the  element  values  in  Figure 
5a  by  G  =  l/Rg,  B  =  -l/coLg,  R  =  Rj,  and  X  =  wL  . 

For  the  open-circuit  tests,  the  low-voltage  winding  was  excited  with  the  high-voltage 
winding  open,  as  shown  in  Figure  5b.  Since  only  limited  current  could  be  supplied  by  the 
480-V  3-phase  (277  V/phase)  commercial  power  hook-up  at  the  test  facility,  the  short-circuit 
tests  for  each  of  the  transformers  were  conducted  by  driving  the  high-voltage  winding  with 
the  low-voltage  winding  shorted.  This  was  done  to  insure  that  the  transformers  were  being 
excited  with  normal  current  levels.  This  change  in  the  test  procedure  requires  a  slight 
modification  of  the  expressions  used  for  computing  the  R  and  X  elements  in  Figure  5c:  the 
turns  ratio  "a"  in  the  expressions  in  Figure  5b  should  be  taken  to  be  unity.  This  modification 
is  needed  because  the  measurement  is  performed  on  the  same  winding  on  the  transformer 
where  the  R  and  X  elements  are  located.  Table  2  presents  the  measured  circuit  values  for  the 
three  transformers. 
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a.  A  Simplified  Transformer  Model. 
Figure  5.  Transformer  Model  and  Test  Set-ups. 
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b.  Test  Configuration  for  Core  Loss  Measurements  (Open  Circuit) 
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Figure  5.  Transformer  Model  and  Test  Set-ups  (Concluded). 
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Table  2 


Measured  Transformer  Circuit  Parameters 


300-kVA  Y-Y  Transformer  T1  (a  =  26) 


Phase 

Rc(Q) 

Lc(h) 

Rtffl) 

(mh) 

1 

157 

0.40 

3.6 

26.3 

2 

273 

0.58 

3.7 

25.0 

3 

317 

0.37 

3.9 

25.4 

75-kVA  Y-Y  Transformer  T2  (a  =  60) 


Phase 

Rc(Q) 

Lc(h) 

Rt(Q) 

Lt  (mh)  1 

1 

53 

0.09 

20.9 

51.5 

2 

136 

0.  14 

20.4 

52.3 

_ 3 _ 

337 

0.09 

21.5 

51.5 

75-kVA  Y-  Transformer  T3  (a  =  60) 


Phase 

Rc(Q) 

Lc(h) 

Rt(«) 

1^  (mh) 

1 

145 

0.70 

22.4 

84.5 

2 

145 

0.81 

19.9 

102.5 

3 

156 

0.87 

23.4 

90.7 

T 


As  discussed  in  ref.[8],  alternate  transformer  representations  are  possible.  Figure  6 
illustrates  one  possible  model  in  which  the  elements  Rt  and  Lj  are  split,  with  part  being 
located  in  the  primary  circuit  and  part  being  put  in  the  secondary  circuit.  Of  course,  the 
fraction  being  put  in  the  primary  circuit  must  be  scaled  by  a  factor  l/a^  , 


9^ 


+ 


S  "2:v 


Secondary 


Figure  6.  An  Alternate,  More  Precise  Transformer  Model. 


The  equivalent  circuit  model  in  Figure  5a  and  the  values  in  Table  2  are  based  on  the 
assumption  of  a  linear  operation  of  the  transformer.  The  nonlinear  core  saturation 
characteristics  of  the  transformers  were  not  measured;  consequently,  the  above  transformer 
model  is  not  sufficient  for  a  complete  nonlinear  characterization  of  the  devices.  As  will  be 
noted  in  Section  7.0,  these  transformers  were  frequently  driven  into  nonlinear  operating 
regions  by  the  dc  neutral  excitation  ctirrents.  It  is  possible  that  with  the  data  presented  in  this 
report,  along  with  a  reasonable  assumption  of  the  transformer  core  behavior,  a  suitable 
nonlinear  operating  model  could  be  developed.  Such  a  model,  however,  is  beyond  the  scope 
of  the  present  effort,  and  is  not  considered  further  in  this  report. 

For  each  of  the  transformer  configurations  depicted  in  Figures  2  to  4,  the  dc  injection 
circuitry  is  shown  in  the  lower  portion  of  the  circuit.  With  the  switch  position  such  that  the 
battery  is  in  the  circuit,  a  dc  current  is  injected  into  the  neutral  circuit.  A  measurement  of  the 
approximate  ac  neutral  current  with  the  battery  in  the  circuit,  but  without  dc  injection,  can  be 
obtained  with  the  switch  in  the  normally  closed  position.  For  dc  injection,  the  switch  is 
moved  to  the  normally  open  position. 
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The  overall  measurement  program  was  divided  into  22  different  tests,  which  are 
summarized  in  Table  3.  This  table  shows  the  test  date,  the  transformer  configuration,  and  the 
injected  dc  current  levels  (I^c).  as  well  as  the  data  file  identification  scheme  used  for 
archiving  the  measured  data.  In  each  test,  measurements  of  both  voltage  and  current  were 
made  at  locations  to  be  discussed  in  detail  below.  As  indicated  by  the  footnote  in  Table  3,  the 
symbol  x  in  the  file  identifier  string  represents  either  I  or  V,  depending  on  the  type  of 
measurement  made.  The  yy  and  zz  symbols  in  the  file  identifier  represent  the  individual  test 
numbers,  01  through  12. 
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Table  3. 

Test  Parameters 


Test 

No. 

File  ID* *• 

Date 

— 

Config. 

0 

lA000xl3.0yy 

6/13/90 

0 

lA 

1 

lAD06xl3.0 

6/13/90 

0.006 

lA 

2 

6/13/90 

0.1 

lA 

3 

lAD60xl3.0 

6/13/90 

0.6 

lA 

4 

■TO^«giki.ga 

6/13/90 

1.0 

lA 

5 

■EQsmEESl 

6/13/90 

1.5 

lA 

6 

lB10Dxl3.0yy 

6/13/90 

10 

IB 

7 

lB25Dxl3.0 

6/13/90 

25 

IB 

8 

6/14/90 

25 

IB 

9 

6/14/90 

3.0 

lA 

10 

6/14/90 

4.5 

lA 

11 

lA5D5xl4.0 

6/14/90 

5.5 

lA 

12 

lA000xl4.Pyy 

6/14/90 

0 

lA 

13 

lA5D5xl4.Pyy 

6/14/90 

5.5 

lA 

14 

2B000xl4.oyy 

6/14/90 

0 

2 

15 

2B5D0xl4.0yy 

6/14/90 

5.0 

2 

16 

2B10Dxl4.0yy 

6/14/90 

10 

2 

17 

2B30Dxl4.0yy 

6/14/90 

30 

2 

18 

2B20Dxl4.0yy 

6/14/90 

20 

2 

19 

2B57Dxl4.0yy 

6/14/90 

57 

2 

21 

2B57Fxl4.0 

YY 

6/14/90 

57 

2 

22 

2B57Fxl5.0 

yy 

16/1 5/90 

57 

2 

Comments 


Base  Case  for  config.  1 A 
26  sec  duration 


26  sec  duration 


26  sec  duration 


Only  six  channels, 
onl 


Incorrect  channel  8  offset 
channel  4  CT  out 


Incorrect  channel  8  offset 
channel  4  CT  out 


Incorrect  channel  8  offset 
channel  4  CT  out 


Incorrect  channel  8  offset 
channel  4  CT  out 


Incorrect  channel  8  offset 
channel  4  CT  out 


Incorrect  channel  8  offset 
channel  4  CT  out 


•  Unless  noted,  all  tests  were  of  15  seconds  duration,  all  12  channels  (Fig.  7)  of  data  were 
taken,  and  the  test  sequence  was  as  follows;  ac  power  on,  data  recorder  on,  dc  injection  on  (2- 
3  seconds),  dc  injection  off  (2-3  seconds),  dc  injection  on  (8-10  seconds),  dc  injection  off, 
data  recorder  off,  ac  power  off.  X  =  I  or  V,  yy  "  01  12,  zz  =  01  06. 

*•  Configuration  is  essentially  that  of  Configuration  IB  with  the  primary  ground  of 
transformer  T2  opened. 
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4.0  INSTRUMENTATION 


The  instrumentation  configuration  for  these  tests  is  shown  in  Figure  7,  and  a  list  of  the 
equipment  used  is  presented  in  Table  4.  The  main  part  of  the  data  acquisition  process 
involved  12  simultaneous  voltage  and  current  measurements  using  the  Astromed  95000 
digital  strip  chart  recorder.  The  line-neutral  voltages  on  the  high-voltage  lines  at  terminals 
A,  B,  and  C  of  transformer  T2  were  measured  on  channels  1,  3,  and  5  of  the  recorder, 
respectively.  The  phase  currents  in  lines  A  and  B  were  measured  in  channels  2  and  4  with  the 
neutral  current  being  monitored  in  channel  6.  The  phase  current  in  line  C  was  not  measured 
because  of  insufficient  data  channels.  On  the  load  (secondary)  side  of  T2,  load  current 
measurements  at  A,  B,  and  C  were  made  on  channels  ^10,  aijd  12,  and  line-neutral  voltage 
measurements  were  made  on  channels  7,  9,  and  1 1 . 

Different  types  of  current  and  voltage  probes  were  used  for  these  measurements,  as 
indicated  in  Figure  7.  These  probes  all  need  a  calibration,  or  "correction,"  factor  to  relate 
their  measured  output  responses  to  the  physical  voltage  or  current  quantity  being  measured. 
These  factors  are  listed  in  Table  4. 

For  most  of  these  probes,  the  probe  conection  factors  are  practically  independent  of 
frequency;  consequently,  a  frequency-independent  probe  correction  factor  is  used.  The  two 
T&M  Research  current  probes  used  on  the  12.47-kV  line  section  for  channels  2  and  4, 
however,  were  found  to  have  a  response  that  varied  appreciably  at  frequencies  below  100  Hz. 
Figme  8  shows  the  measured  transfer  impedance  magnitude  of  these  probes.  This  quantity  is 
defined  as  the  induced  voltage  across  the  probe,  Vp  ,  when  measured  by  a  voltmeter  with  a 
50-0  input  impedance,  Z^,  created  by  a  current  of  amps,  and  is  expressed  as 

Z,  =  |Vp  /  ^ol^.joQ-  curve  in  Figure  8  indicates  that  the  lower  3-dB  point  for  these 

probes  is  approximately  at  70  Hz.  Thus,  for  these  probes  to  be  used  accurately  at  60  Hz,  a 
correction  for  this  frequency  fall-off  must  be  used.  This  problem  is  discussed  in  more  detail 
later  in  the  report.  Note  that  at  frequencies  over  about  200  Hz,  the  probe  transfer  impedance 
is  a  constant  of  about  0.45  0  . 
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Test  Instrumentation. 


Table  4. 
Test  Equipment 


Equipment 

Model  No 

Astromed 

Measurement 

Channels 

Correction  Factor 

Dranetz  Power  Demand 
Analyzer 

808 

N/A 

Dranetz  Current  Probe 

2013 

— 

Times  3  accounted  for  by 
Power  Demand  Analyzer 

Tektronics  Voltage 
Probe 

6015 

1.3,5 

Times  1000 

T&M  Research  Current 
Viewing  Resistor 
(CVR) 

1M-2N 

6 

Times  100 

T«S:M  Research  Current 
Viewing  Resistor 
fCVR) 

F-5000-20N 

8,10,12 

Times  1000 

T&M  Research  Current 
Sensor 

CT6.5.685- 

5-100 

2,4 

Times  1000  with  R  =  0.10 
Times  2000  with  R  =  0.050 
Times  2  with  R  =  50  0 
(Frequency  correction 
needed  for  500  termination) 

Tektronics  Voltage 
Probe 

6062B 

7,9 

Times  10 

Tektronics  Voltage 
Probe 

6105 

11 

Times  10 

Astromed  12  Channel 
Digital  Strip  Chart 
Recorder  with  Deep 
Memory 

95000 

N/A  1 

15 


60  Hi 


Figure  8.  Measured  Transfer  Impedance  (Zj)  for  T&M  Research  Current  Probe. 


The  same  probes  were  used  occasionally  to  monitor  the  input  current  to  the  test  setup, 
but  in  this  case,  optional  low-resistance  terminators  on  the  probes  were  used.  While  this 
reduces  the  sensitivity  of  the  probe  (i.e.,  lowers  the  transfer  impedance),  it  extends  the  lower 
3-dB  point  well  below  60  Hz.  Thus,  the  frequency-independent  scalars  for  these  probes, 
which  are  supplied  in  the  data  records,  are  sufficient  for  the  probe  calibration  factors  in  this 
case. 


For  each  test,  the  input  voltage,  current,  and  power  were  monitored  by  the  Dranetz 
Power  Demand  Analyzer,  which  provided  readings  of  these  quantities  on  a  data  tape.  These 
quantities  served  as  a  check  of  the  overall  operation  of  the  system,  and  provided  an  additional 
means  of  correcting  the  V  and  I  measurements  on  the  high-voltage  portion  of  the  lines  to 
insure  that  a  power  balance  is  achieved.  As  will  be  discussed  later,  these  measmements 
provide  a  way  of  insuring  that  the  probe  corrections  are  done  properly. 
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5.0  TEST  PROCEDURES 


For  each  of  the  transformer  configurations  shown  in  Figures  2  through  4,  the  tests 
indicated  in  Table  3  were  conducted.  Once  the  transfonners  had  been  connected  in  the  proper 
manner,  the  rheostat  in  the  dc  circuit  was  adjusted,  with  no  ac  excitation  on  the  transfoimers, 
to  provide  the  desired  dc  neutral  current  level.  The  battery  switch  was  then  placed  in  the 
normally-closed  (NC)  position  so  no  dc  excitation  of  the  transformer  neutral  was  present,  and 
the  transformers  were  energized  with  the  commercial  480-V  power.  The  ac  neutral  cunent 
was  measured  for  the  case  of  no  dc  excitation.  After  this  was  complete,  the  battery  was 
switched  into  the  circuit  by  energizing  the  12  V  control  relay.  The  complete  set  of  voltage 
and  current  measurements  was  made  for  a  fixed  period  of  time  (13  seconds  for  most 
measurements). 

Since  the  time  required  to  reach  a  steady  state  response  after  the  power  was  switched 
on  was  found  to  be  significantly  less  than  this  measurement  time,  the  15  second  measurement 
time  was  more  than  adequate  for  this  test.  Furthermore,  in  most  of  the  tests,  it  was  possible  to 
disconnect  and  then  reenergize  the  dc  circuit,  thereby  providing  a  clear  view  of  the  transient 
relaxation  time  of  the  response. 

While  damage  to  the  transformer  was  permissible  during  this  test  program,  the 
duration  of  the  dc  current  injection  was  sufficiently  short  to  preclude  such  damage,  even  in 
the  most  severe  case  of  current  injection.  Consequently,  no  damage  to  the  transformers  or 
loads  was  noted,  aside  from  the  blowing  of  several  fuses  on  the  power  feed  line  to  the  test 
setup.  This  occurred  for  the  initial  in-rush  current  of  the  energizing  transfoimers,  but  not  for 
the  dc  current  injection. 
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6.0  DATA  PROCESSDJG  FOR  THE  DISTRIBUTION  TRANSFORMER  TEST 
6.1  Overview 

As  indicated  in  Section  4,  the  voltage  and  current  data  from  this  test  program  were 
measured  on  a  12-channel  Astromed  MT-9500  transient  digitizer.  The  digitizer  provided  a 
plot  of  the  time  histories  of  each  of  the  measured  results  for  a  short  period  of  time  at  the 
beginning  of  each  measurement,  as  well  as  a  binary  output  stream  of  data  which  was  captured 
onto  a  PC  disk.  As  indicated  in  Table  3,  data  were  acquired  for  26  seconds  for  the  first  three 
transformer  configurations  tested  (lAOOO,  1AD06  and  lADlO).  The  raw  data  from  the 
digitizer  was  contained  in  32  packets  of  binary  data,  each  containing  4096  data  samples.  This 
gave  a  total  of  131,072  individual  data  points,  with  a  resulting  raw  data  file  size  of  266 
KBytes  for  each  of  the  12  measurements  for  these  configurations.  The  data  acquisition  for 
the  remaining  19  test  configurations  lasted  for  13  seconds;  consequently,  the  raw  data  file  size 
was  133  KBytes  each.  For  all  cases,  the  resulting  time  interval  between  data  points  was  t  = 
0.2  milliseconds,  or  1.2%  of  the  60-Hz  period.  From  this  series  of  tests,  252  separate 
transient  waveforms  were  measured. 

A  data  translation  program  -.  as  provided  by  MRC  [9]  to  read  each  of  the  raw  binary 
files  and  convert  them  to  ASCII  files  of  ordered  time  and  response  data  pairs.  The  resulting 
sizes  of  the  ASCII  data  files  were  3.54  MB>tes  for  me  data  taken  over  the  26  sec  time  period 
and  1.8  MBytes  for  the  13  sec  data.  This  resulted  in  approximately  538  MBytes  of  ASCII 
data  for  this  test. 

The  first  step  in  the  analysis  procedure  for  this  test  was  to  perform  a  "quick-look" 
(overview)  of  the  measured  data.  This  was  done  to  see  if  there  was  any  evidence  of 
transformer  saturation  for  the  different  levels  of  dc  injection.  This  processing  was  done  for 
all  of  the  data  records.  After  this  preliminary  processing  was  completed  and  the  results 
analyzed,  a  more  detailed  spectral  analysis  of  the  data  was  undertaken.  This  analysis  resulted 
in  plots  of  the  time  variations  of  the  60-Hz  components  of  the  measured  responses  along  with 
the  next  5  higher  harmonics.  A  final  analysis  was  conducted  to  determine  the  time  variations 
of  the  real  and  imaginary  components  of  the  60-Hz  power  flowing  along  the  distribution  line. 
This  is  important  in  considering  possible  increased  reactive  power  demand  on  the  generation 
facilities,  because  of  the  MHD-EMP  excitation  of  transformers. 

Because  of  the  large  amount  of  data  acquired  in  this  test  and  the  fact  that  the 
processing  was  done  on  an  IBM  PC-class  machine,  care  was  used  to  develop  these  data 
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processing  routines  in  such  a  way  that  they  would  not  take  an  exorbitant  amount  of 
processing  time  or  use  more  disk  or  memory  space  than  absolutely  required.  In  the  following 
sections,  each  of  these  analyses  is  discussed  in  more  detail,  with  examples  of  the  processed 
data  being  shown. 

6.2  "QUICK-LOOK"  ANALYSIS 

Using  the  MRC-provided  data  conversion  program,  a  separate  routine  was  written  to 
read  the  raw  binary  data  files,  convert  them  into  ASCII,  and  generate  a  plot  file  showing  the 
envelope  of  the  measured  voltage  or  current  responses.  For  this  quick-look  analysis,  only  the 
appropriate  scalar  (frequency  independent)  measurement  probe  transfer  functions  were  used. 
These  parameters  were  imbedded  in  the  binary  data  file  and  were  used  to  obtain  the  measured 
responses  in  engineering  units  (i.e.,  in  volts  and  amps).  The  requirement  for  an 
implementation  of  frequency-dependent  probe  corrections  is  discussed  in  the  following 
section  on  harmonics.  Because  of  the  large  storage  requirements  for  retaining  the  converted 
ASCII  data,  these  data  files  were  discarded  once  they  were  plotted. 

As  an  example  of  results  from  this  phase  of  the  analysis.  Figure  9  illustrates  typical 
output  data  for  measured  currents  on  channel  2  (the  high-voltage  side)  and  channel  8  (the 
low-voltage  side)  of  the  first  transformer.  These  are  the  actual  transient  waveforms,  shown 
over  a  limited  time  window,  for  the  case  of  no  dc  injection  into  the  transformer  (Fig.  9a,  for 
Configuration  lAOOO)  and  the  case  of  5.5  A  injection  (Fig.9b  for  Configuration  1A5D5). 
Part  a  of  the  figure  is  for  no  dc  current  injection,  and  exhibits  a  reasonably  clean  60-Hz 
sinusoidal  waveform.  Part  b  on  the  other  hand,  illustrates  a  highly  perturbed  waveform 
arising  from  the  saturation  of  the  transformer  cores  due  to  the  dc  current  injection.  The 
secondary  current  for  no  dc  injections  is  shown  in  part  c.  Note  that  some  distortion  is  present 
on  the  peaks  of  the  normal  60-Hz  sinusoidal  waveform.  This  may  indicate  that  the 
transformer  was  near  magnetic  saturation  under  normal  operating  conditions.  The  waveform 
in  part  d  is  clearly  more  distorted  than  the  corresponding  primary  waveform  in  part  b  with 
5.5  A  dc  injection.  The  transformer  shows  no  tendency  to  "block"  these  lower-order 
harmonics  and  may  have  enhanced  them. 
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c. 


0.60  0.62  0.64  0.66  0.68  0.70 

Time  (Sec) 


0.60  0.62  0.64  0.66  0.68  0.70 

Time  (Sec) 


a.  Primary  Current  (Channel  2)  for  0.0  A  Dc  b.  Primary  Current  (Channel  2)  for  5.5  A  Dc 


0.60  0.62  0.64  0.66  0.68  0.70 

Time  (Sec) 


0.60  0.62  0.64  0.66  0.68  0.70 

Time  (Sec) 


Secondary  Current  (Channel  8)  for  0.0  A  Dc  d.  Secondary  Current  (Channel  8)  for  5.5  A  Dc 


Figure  9.  Examples  of  Measured  Transient  Currents  for  Different  Dc  Levels. 
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The  quick-look  analysis  consisted  of  plotting  all  the  da^  files  in  such  a  way  to  show 
the  behavior  of  the  envelope  of  the  curves.  In  this  manner,  the  overall  shape  of  the  waveform 
amplitudes  can  be  determined  and  the  effects  of  saturation  can  be  noted  as  a  fimction  of  the 
dc  current  levels.  As  an  example  of  these  results.  Figure  10  illustrates  the  envelope  curves  for 
the  phase  conductor  current  of  Figure  9b.  Although  the  details  of  the  current  waveshape  are 
lost  in  these  envelope  plots,  the  initial  tum-on  of  the  dc  excitation  is  clearly  evident,  with  the 
ciurent  peaks  growing  to  about  5.2  A  with  an  associated  time  constant  of  the  system  on  the 
order  of  one  second.  Later,  at  about  t  =  2  seconds,  the  dc  is  switched  off,  and  the  envelope 
curve  is  noted  to  drop  back  to  about  2.0  A,  a  value  consistent  with  the  unexcited  case  of 
Figure  9a.  At  about  3.5  seconds  the  dc  injection  is  again  turned  on,  and  the  response  again 
grows  and  becomes  nonlinear. 


0  5  10  15 

Time  (Sec) 


Figure  10.  Example  of  the  Envelope  Plots  for  the  Preliminary  Analysis. 
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6.3  Harmonic  Analysis 


The  idea  behind  performing  a  harmonic  content  analysis  of  the  measured  data  is  to 
understand  the  time-dependent  development  of  the  60-Hz  harmonics  as  the  dc  current 
injection  is  applied  to  the  system.  Clearly,  the  harmonic  content  of  the  waveform  in  Figure 
9a  will  be  significantly  lower  than  in  Figure  9b,  and  as  the  transformer  transitions  from  one 
state  to  another,  the  harmonic  content  must  undergo  a  smooth  transition.  It  is  useful  to  plot 
the  transient  behavior  of  these  harmonics  and  show  the  rate  of  saturation  of  the  transformer. 

Figure  11a  shows  the  Fourier  spectrum  magnitude  of  the  time  domain  current  of 
Figure  9a.  This  was  computed  by  using  a  discrete  fast  Fourier  transform  (FFT)  involving 
2048  time  samples.  It  is  presented  as  a  continuous  spectral  density,  with  a  vertical  scale  in 
units  of  A/Hz.  Normally,  a  periodic  function  of  infinite  duration  with  a  fundamental 
frequency  of  60  Hz  like  those  in  Figure  9  will  have  a  set  of  discrete  spectral  lines  whose 
amplitudes  are  the  coefficients  of  a  Fourier  series.  Because  the  FFT  calculation  uses  a 
windowed  time  function,  the  resulting  FFT  spectrum  consists  of  a  series  of  resonances  of 
finite  width.  As  the  length  of  the  time  record  increases,  the  width  of  these  spectral  lines 
decreases  until  true  impulses  are  obtained  for  a  record  of  infinite  length. 

To  relate  these  calculated  spectra  to  a  Fourier  series  equivalent,  it  is  necessary  to 
multiply  them  by  a  factor  "f,  which  is  the  sampling  interval  in  the  frequency  domain  [10]. 
This  factor  is  inversely  related  to  the  maximum  time  of  the  transient  waveform,  and  is  equal 
to  2.44  Hz  in  the  present  case.  Therefore,  at  the  fundamental  resonant  frequency  of  Figure 
11a,  the  spectral  component  is  about  0.9  A.  If  the  waveform  of  Figure  9a  were  analyzed 
analytically,  assuming  a  perfect  sinusoidal  waveform,  its  60-Hz  Fourier  series  coefficient 
would  be  about  1.0.  This  process  serves  as  a  check  that  the  nmnerically  calculated  Foiuier 
spectrum  is  being  done  properly. 

The  corresponding  spectral  density  magnitude  for  the  current  of  Figure  9b  is  shown  in 
Figure  1  lb.  Here  it  is  clear  that  the  nonlinear  behavior  of  the  transformer  has  resulted  in  a 
waveform  rich  in  spectral  harmonics,  occurring  at  multiples  of  60  Hz. 

To  process  the  entire  13  seconds  of  transient  data,  a  sliding  FFT  window  method  was 
used.  The  FFT  window  consisted  of  1024  sample  points,  with  a  total  period  of  0.204 
seconds.  This  window  started  at  the  beginning  of  the  waveform,  i.e.,  at  t  =  0  ,  and  moved 
along  until  the  end  of  the  total  waveform  was  reached.  This  resulted  in  a  total  of  64  FFTs 
performed  for  each  waveform.  At  each  time  step  in  the  waveform,  the  Fourier  series  spectral 
lines  at  60,  120,  180, 240,  300,  and  360  Hz  were  computed  and  saved. 
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Figure  11.  Examples  of  Computed  Spectra  for  Measured  Responses 
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The  quick-look  data  processing  did  not  include  frequency-dependent  probe 
corrections  for  to  the  low-frequency  fall-off  of  the  probe  transfer  impedances  shown  in  Figure 
8.  Most  of  the  voltage  or  current  probes  used  in  these  tests  were  described  by  frequency- 
independent  calibration  factors  which  simply  scaled  the  measured  probe  responses.  These 
factors  were  incorporated  in  the  generation  of  the  ASCII  data  files  and  subsequently  used  in 
the  data  processing.  For  data  channels  2  and  4,  however,  the  T&M  Research  current  probes 
were  used,  and,  as  indicated  previously  in  Figure  8,  corrections  to  the  probe  transfer 
impedance  were  required  at  frequencies  below  200  Hz  if  an  accurate  determination  of  the 
harmonic  content  of  the  responses  was  to  be  obtained. 

One  way  of  correcting  the  probe  transfer  impedance  is  to  simply  digitize  the 
impedance  shown  in  the  figure  and  provide  a  calibration  table  of  probe  impedance  vs. 
frequency.  Given  a  particular  measured  probe  response,  which  is  to  be  corrected  at  a 
frequency  fg  ,  the  appropriate  complex-valued  probe  correction  factor  is  interpolated  from 
the  data  table  and  then  divided  into  the  probe  response  to  obtain  the  corrected  current.  Note 
that  this  correction  factor  must  involve  both  the  magnitude  and  the  phase  of  the  probe  transfer 
impedance,  although  only  the  magnitude  function  is  shown  in  the  figure.  This  approach  is 
identical  to  many  of  the  instrumentation  probe  correction  procedures  used  for  fast  transient 
measurements  for  HEMP  testing. 

An  alternative  to  interpolating  from  the  data  table  is  to  use  an  analytical  curve-fit  to 
the  measured  probe  transfer  impedance  function.  An  examination  of  the  spectral  magnitude 
function  for  the  probe  in  Figure  8  indicated  that  the  probe  is  behaving  as  a  1-pole  high-pass 
filter  at  low  frequencies.  Using  a  network  analyzer,  a  precise  measurement  of  the  lower  3-dB 
frequency  of  the  probe  was  made  by  MRC  using  a  variable  frequency  marker,  giving  a  value 
of  116.8  Hz,  or  about  740  radians/second.  Using  this  low-frequency  break-point  and  the 
knowledge  that  the  high-frequency  asymptote  of  the  probe  impedance  is  about  0.45  0  ,  the 
following  1-pole  fit  can  be  developed  for  the  probe  transfer  impedance: 


Z.(Q)) 


(ja)x0.45) 
(jo)  +  740.) 


(fl) 


(1) 


where  u  is  the  angular  frequency,  <a  =  2Tf.  In  the  actual  probe,  there  is  an  additional  high- 
frequency  fall-off,  along  with  a  possible  probe  resonance  at  10  MHz,  but  because  we  are 
interested  only  in  the  low-frequency  behavior  of  the  probe,  these  high-frequency  effects  are 
not  modeled  here. 
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The  probe  correction  factor  defined  in  Eq.  (1)  above  is  a  complex-valued  function, 
having  both  a  magnitude  and  a  phase  compunent.  Figures  12a  and  12b  present  these 
quantities  as  determined  from  the  above  expression,  both  shown  as  a  function  of  frequency. 
The  magnitude  of  at  f  =  60  Hz  is  seen  to  be  about  0.20  0  ,  which  confirms  the  validity  of 
Eq.  (1). 


a.  Magnitude 
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b.  Phase  Angle 

Figure  12.  Plots  of  the  Analytical  Current  Probe  Impedance  Fimction. 
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Figure  8  shows  that  the  magnitude  of  at  the  same  frequency  is  slightly  different, 
having  a  measured  value  of  about  0.26  0.  This  difference  arises  from  the  fact  that  the  3-dB 
frequency  point  of  f3jjB  »  90  Hz  taken  from  the  curve  in  Figure  8  is  different  from  that  of 
fjdB  •  1  16.8  Hz  obtained  by  measuring  this  frequency  directly  with  the  marker  in  the 
network  analyzer.  The  discrepancies  between  these  two  3-dB  frequencies  indicate  that  there 
is  also  some  uncertainty  in  the  probe  calibration  process,  a  fact  that  will  be  evident  when  an 
energy  balance  is  performed  later  in  Section  6.4  of  this  report. 

As  an  example  of  the  effect  of  the  frequency  dependence  of  this  probe  correction 
function,  the  response  shown  in  Figure  1  lb  has  been  reprocessed  using  the  probe  transfer 
impedance  of  Eq.  (1).  Figure  13  shows  the  resulting  current  spectrum,  corrected  for  the  low- 
frequency  probe  fall-off.  This  spectrum  exhibits  a  definite  increase  in  the  low-frequency 
response,  as  expected  from  the  probe  characteristics. 

To  illustrate  this  probe  correction  procedure  on  the  measured  spectral  harmonics  of 
the  line  current,  the  time  history  of  the  primary  60-Hz  spectral  coefficient  and  the  next  5 
harmonics  is  shown  in  Figure  14a  for  the  case  of  no  dc  injection  into  the  transformer.  The 
primary-side  current  is  significantly  larger  at  60  Hz  than  at  the  higher  harmonics,  and  the 
relative  ratios  of  the  harmonics  appear  to  be  constant  as  time  progresses.  For  the  case  with 
5.5  A  of  dc  injection,  Figure  14b  shows  the  associated  spectral  content.  As  discussed  earlier, 
the  dc  injection  was  switched  off  at  about  t  =  2  seconds,  and  turned  on  again  at  about  t  =  3.5 
seconds.  The  effect  of  this  procedure  on  the  spectral  content  is  clearly  evident.  Note  that  the 
60-Hz  spectral  component  is  slightly  affected  by  this  dc  excitation  (see  the  slight  dip  at  t  =  2 
seconds). 

6.4  Power  Analysis 

As  mentioned  previously,  simultaneous  voltage  and  current  measurements  were  made 
on  phases  1  and  2  of  the  high-voltage  side  of  the  first  transformer  and  on  all  three  phases  of 
the  terminating  resistors  of  the  loads.  Consequently,  these  measurements  can  be  used  to  infer 
the  amount  of  complex  power  flowing  on  the  lines  and  into  the  load.  Unfortunately,  for  the 
high-voltage  section  on  phase  3,  only  the  voltage  was  measured,  so  its  power  cannot  be 
determined  directly.  In  addition  to  making  these  power  measurements,  we  placed  a  Dranetz 
Power  Demand  Analyzer  at  the  source  end  of  the  transformers,  and  were  thereby  able  to 
monitor  both  real  and  reactive  powers  on  all  phases. 
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Figure  13.  Corrected  Spectrum  of  Fig.  9b  Using  the  Probe  Correction  Function  of  Fig.  12. 
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Figure  14.  Examples  of  Principal  Harmonic  Content  for  Measured  Waveforms. 


At  a  particular  frequency  w  the  complex  power  is  defined  in  terms  of  the  peak  voltage 
and  current  phasor  components  as 

P(w)  -  1/2 1(o))  V(a))  (W) .  (2) 

Relating  these  phasors  to  the  computed  spectral  component  at  60  Hz  permits  the 
estimation  of  the  delivered  power  at  60  Hz.  Note  that  Eq.  (2)  defines  a  complex  power,  the 
value  of  which  is  critically  dependent  on  the  phase  relationship  between  the  measured  voltage 
and  the  current. 

Applying  Eq.  (2)  to  the  voltage  and  current  measurements  at  the  208-V  side  of  the  T2 
transformer  provided  estimates  of  the  power  delivered  to  the  loads.  These  could  be  compared 
with  the  input  power  measurements  from  the  power  demand  meter  in  the  source  circuit. 
Good  agreement  in  these  comparisons  was  noted. 

However,  for  the  checks  on  the  power  flow  on  phases  1  and  2  in  the  high-voltage 
section  of  the  line,  between  T1  and  T2,  it  was  found  that  the  agreement  was  not  satisfactory. 
As  an  example,  Figure  15a  presents  the  measured  60  Hz  power,  both  real  (watts)  and 
imaginary  (vars),  as  a  function  of  time  for  the  0.0  A  dc  injection  case  for  configuration  lA. 
The  power  demand  meter  in  this  case  registered  a  real  power  of  approximately  19.1  kW,  with 
0.0  var.  The  real  power  in  this  figure  is  approximately  correct,  but  the  reactive  power  is  too 
large.  This  error  arises  due  to  the  uncertainty  in  the  measured  data  for  the  low-frequency 
probe  correction,  as  discussed  previously. 

In  the  absence  of  additional  probe  characterization  data,  a  remedy  to  this  problem  was 
sought  by  modifying  the  probe  transfer  impedance  so  as  to  provide  the  proper  power 
calculations  at  60  Hz.  This  was  done  by  multiplying  the  probe  transfer  impedance  by  the 
complex  quantity  (0.90  -  j  0.37).  Figure  15b  illustrates  the  resulting  60  Hz  power  in  this  case. 
Note  that  the  computed  reactive  power  is  now  consistent  with  that  measured  from  the  power 
demand  meter. 

Figure  16  presents  the  computed  60-Hz  power  as  a  function  of  time  for  phase  1 
(channels  1  and  2)  for  the  case  of  the  5.5-A  dc  injection  that  has  been  discussed  previously. 
Notice  that  the  presence  of  the  dc  excitation  places  a  very  significant  reactive  power  demand 
on  the  generator  feeding  the  transformers. 
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Figure  15.  Complex  60-Hz  Power  Computed  for  Phase  1  with  No  Dc  Excitation. 
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Figure  16.  Complex  Power  for  Phase  1  Configuration  lA,  With  5.5  A  Dc  Excitation 
Switched  On,  Off,  and  On. 
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7.0  OBSERVATIONS  FROM  THE  TEST  DATA 


From  the  data  taken  in  this  test,  a  number  of  observations  regarding  transformer  op¬ 
eration  under  a  combined  dc  and  ac  excitation  can  be  made.  This  section  will  summarize 
some  of  the  important  observations  for  the  data  and  will  discuss  selected  highlights  of  the 
responses. 

7.1  Configuration  lA 

The  first  case  to  be  discussed  is  configuration  1  A.  As  shown  in  Figure  2,  this  consists 
of  two  Y-Y  transformers  with  the  neutrals  grounded,  and  with  the  possibility  that  a  dc  current 
could  be  injected  in  the  neutral  of  the  primary  circuit.  It  is  of  interest  to  observe  the  behavior 
of  the  phase  current  and  the  phase-neutral  voltage  on  the  12.47-kV  section  of  the  line  for 
varying  levels  of  dc  injection.  For  no  dc  injection,  as  noted  in  Figure  9,  the  line  current 
response  consists  simply  of  the  normal  6O-H2  phase  current.  As  the  dc  current  level  increases, 
the  line  current  exhibits  abnormal  behavior.  Figure  17  shows  the  envelope  plots  for  the  line 
current  on  phase  1  for  dc  current  injections  of  0.0,  0.6,  1.0,  1.5,  3.0,  and  5.5  A.  As  discussed 
in  the  previous  section,  these  envelope  curves  have  been  determined  using  the  scalar  (i.e., 
frequency  independent)  current  probe  calibration  factors.  As  such,  they  can  only  provide  a 
qualitative  indication  of  the  current  behavior. 

For  these  injection  tests,  the  ac  power  to  the  transformer  was  turned  on,  and  several 
seconds  later  the  dc  battery  source  was  applied.  The  data  recorder  was  activated  after  the  dc 
was  applied,  so  that  the  complete  tum-on  transient  response  due  to  the  dc  injection  is  not 
always  evident  in  the  measurements.  Later,  two  to  three  seconds  into  the  test,  the  dc  power 
was  removed  from  the  neutral  circuit  in  such  a  way  that  neutral  current  could  continue  to 
flow.  This  dc  excitation  was  then  restored  several  seconds  later. 

The  5.5-A  case  in  Figure  17  shows  clearly  the  various  stages  of  the  response  caused 
by  applying  the  dc  excitation.  The  initial  tum-on  transient  is  evident,  as  is  the  tum-off  of  the 
excitation  at  about  t  =  2  seconds.  It  is  interesting  to  note  that  there  appears  to  be  an  L/R  time 
constant  associated  with  this  overall  envelope  shape  which  is  a  function  of  the  dc  current 
injection  level.  The  observed  relaxation  times  of  this  envelope  appear  to  be  much  shorter  as 
the  dc  levels  increase.  With  the  ac  power  removed,  as  seen  in  the  1.0-A  case  where  an  ac  fuse 
blew,  the  system  response  dies  very  quickly  to  zero  with  no  ac  excitation.  In  this  case,  the 
expected  1/3  of  the  1.0-A  dc  injection  current  flowing  in  the  phase  1  conductor  is  not 
’■egistered,  because  of  the  low-frequency  fall-off  of  the  probe  response  and  its  subsequent 
..lability  to  accurately  measure  dc  currents. 
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Figure  17.  Envelope  Plots  of  the  12.47-kV  Phase  1  Current  for  Configuration  lA  for 
Different  Levels  of  Dc  Injection  Switched  On,  Off,  and  On. 
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Figure  18  presents  the  corresponding  phase  1  to  neutral  voltage  (Channel  1)  on  the 
12.47-kV  section  of  line  for  the  six  different  levels  of  dc  current  injection.  There  is  a 
noticeable  change  of  this  quantity  as  a  function  of  the  dc  levels,  but  this  change  is  not  as  large 
as  the  change  in  the  line  current. 

The  3.0-A  case  shown  in  Figure  17  exhibits  a  noticeable  asymmetry  about  the  zero 
cxiirent  axis.  As  the  dc  level  increases,  the  degree  of  asymmetry  decreases  but  is  still  present. 
This  suggests  the  possibility  that  each  phase  of  the  line  has  a  different  response.  Figure  19 
shows  a  comparison  of  the  phase  1  and  phase  2  currents  (measured  in  Channels  2  and  4, 
respectively)  for  0.0,  3.0  and  5.5  A  dc  injection.  Unfortunately,  the  phase  3  currents  were  not 
measured.  As  noted  in  the  zero  dc  excitation  case,  there  is  a  considerable  imbalance  of  the  ac 
phase  current  envelopes  with  a  noticeable  baseline  shift  for  the  phase  2  response.  For  the  3.0- 
A  injection,  there  is  a  distinct  difference  between  the  two-phase  current  envelopes,  but  as  the 
dc  current  level  increases,  these  envelopes  become  more  alike.  In  theory,  this  effect  might  be 
caused  by  the  different  saturation  states  of  the  three  transformer  cores  for  low  levels  of 
applied  dc  excitation. 

The  behavior  of  the  neutral  current  envelopes  from  Channel  6  is  shown  in  Figure  20 
for  the  six  levels  of  dc  injection.  It  is  clear  for  under  no  dc  excitation  the  neutral  current  is 
quite  small,  indicating  that  the  three  phases  of  the  power  system  are  well  balanced.  With  the 
application  of  a  small  amount  of  dc  however,  there  is  a  marked  increase  in  neutral  current. 
Even  a  0.6-A  dc  injection  causes  peak  neutral  currents  on  the  order  of  8-A  to  flow.  This 
increase  in  the  neutral  current  is  presumably  caused  by  the  different  saturation  of  the 
transformer  cores,  which  causes  a  system  imbalance.  In  addition,  the  harmonic  content  for 
these  neutral  currents  is  significantly  higher  than  that  for  the  currents  on  the  phases.  In  fact, 
with  the  exception  of  the  0.0- A  dc  case,  the  harmonics  in  the  neutral  current  at  180  and  300 
Hz  were  larger  than  the  fundamental  component  at  60  Hz.  It  is  interesting  to  note  that  in 
these  cases,  the  peak  value  of  the  neutral  current  is  significantly  larger  than  either  the  peaks  of 
the  phase  conductor  currents  or  the  injected  dc  current.  This  larger  peak  is  caused  by  the  high 
harmonic  content  in  the  neutral  current  waveform,  which  arises  from  the  core  saturation. 
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Figure  19.  Envelope  Plots  for  Phases  1  and  2  Currents  for  Various  Levels  of  Dc  Injection 
Switched  On,  Off,  and  On  (Phase  3  Not  Measured). 
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Figure  20.  Envelope  Plots  for  the  Neutral  Currents  in  the  12.47-kV  Line  for  Various  Levels 
of  Dc  Injection  Switched  On,  Off,  and  On. 


Figure  21  illustrates  the  behavior  of  the  load  current  envelopes  at  the  secondary  side 
of  the  second  transformer  for  0.0  and  5.5  A  dc  injection.  Although  the  load  currents  are  more 
or  less  balanced  when  there  is  no  dc  injected  into  the  transformers,  the  load  currents  on  the 
208  V  side  become  unbalanced  in  the  presence  of  the  dc  excitation.  This  observation,  as  well 
as  that  concerning  the  neutral  current,  suggests  that  the  saturation  characteristics  of  each 
phase  of  the  transformer  core  are  slightly  different. 

Spectral  Analysis 

As  mentioned  in  the  previous  section,  the  behavior  of  the  envelope  curves  does  not 
tell  the  complete  story  of  the  voltage  or  current  behavior,  because  these  waveforms  consist  of 
a  number  of  harmonic  components  as  well  as  the  60-Hz  fundamental  component.  Using  the 
spectral  processing  method  discussed  previously,  the  harmonic  content  of  the  12.47-kV  phase 
1  current  and  voltage  has  been  computed  and  is  plotted  in  Figure  22  as  a  function  of  the  dc 
cuaent  injection.  As  was  noted  in  Figure  18,  the  harmonic  content  of  the  line  voltage  is  not 
as  sensitive  to  the  dc  current  injection  as  is  the  line  current,  which  exhibits  a  continual 
increase  in  all  of  the  harmonics  as  the  dc  injection  increases.  The  60-Hz  fundamental 
component  of  both  the  current  and  voltage  appears  to  be  rather  insensitive  to  the  dc  injection, 
however. 

Reactive  Power  Demand 

It  is  of  interest  to  determine  the  change  in  reactive  power  flowing  in  the  system,  and 
the  especially  the  possible  increase  in  the  reactive  demand  placed  on  the  generation  facilities 
by  an  increase  in  the  dc  injection  level.  This  amounts  to  determining  the  change  in  the  phase 
angles  between  the  voltage  and  current  60-Hz  spectral  components  as  the  dc  level  varies. 
From  the  complex  power  computed  using  the  phase  1  voltage  and  current  spectra,  the 
increase  in  reactive  power  on  phase  1  is  illustrated  in  Figure  23a.  There  is  a  nearly  linear  fit 
to  the  reactive  demand  as  a  function  of  the  dc  injection.  At  the  injection  level  of  5.5  A,  the 
reactive  power  is  about  15  kvar.  For  the  0.0  A  dc  injection,  the  normal  power  per  phase  is 
about  18  kW.  This  increase  in  reactive  power  can  pose  a  considerable  problem  for  come 
generation  equipment. 

Figure  23b  shows  the  behavior  of  the  reactive  power  for  very  small  levels  of  dc 
injection.  For  the  cases  of  0.0,  0.06  and  0.1  A  of  dc  injection,  the  reactive  power  increase 
was  essentially  zero.  Above  0.1  A,  the  reactive  power  begins  to  increase  linearly.  Thus,  it 
appears  that  the  limit  of  linear  operation  for  this  particular  configuration  is  at  about  0.1  A  dc. 
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Figure  2 1 .  Load  Current  Envelopes  for  Configuration  1 A  with  Dc  Switched  On.  Off.  and  On 
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Figure  22.  Harmonic  Content  vs.  Dc  Current  Level  for  Phase  1,  Configuration  lA. 
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Figure  23.  Reactive  Power  on  Phase  1  vs.  Dc  Current  Injection  Level. 
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As  noted  in  the  envelope  plots,  there  appears  to  be  an  L/R  time  constant  associated 
with  the  system,  and  it  appears  that  this  quantity  is  a  function  of  the  dc  current  injection  level. 
This  time  constant  could  be  determined  directly  from  the  envelope  plots,  or  from  the  transient 
behavior  of  the  reactive  power  demand,  as  the  dc  injection  is  switched  off.  This  latter 
approach  was  used:  it  involved  finding  the  time  taken  for  the  reactive  power  to  drop  to  0.368 
of  its  steady-state  value  with  the  dc  injection  on.  This  corresponds  to  a  1/e  time  constant,  a 
commonly  used  measure  in  linear  systems.  In  the  present  case,  with  nonlinear  core 
saturation,  it  is  by  no  means  certain  that  an  exponential  decay  of  these  responses  is  indicated. 
Nevertheless,  this  provides  a  convenient  way  of  describing  the  transient  decay  of  the  system. 

Figure  24  presents  this  observed  time  constant  for  the  1 A  configuration  and  illustrates 
a  decreasing  system  relaxation  time  as  the  dc  injection  increases.  This  finding  is  consistent 
with  the  idea  that  the  normally  large  core  inductance  decreases  drastically  as  the  core 
saturates,  thereby  decreasing  the  effective  L/R  time  constant  of  the  system.  In  the  case  of  no 
dc  injection,  there  normally  will  be  an  in-rush  transient  which  is  observed  when  the 
transformer  is  first  energized  by  the  ac  power.  Such  transients  can  last  for  seconds.  No  such 
transient  was  observed  in  the  data  for  this  configuration  with  zero  dc  injection. 


Figure  24.  System  Relaxation  Time  vs.  Dc  Current  Injection  Level  for  Configuration  1  A. 
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7.2  Conflguration  2 


Configuration  2,  as  shown  in  Figure  4,  consists  of  a  Y-Y  transformer  at  the  source  end 
of  the  12.47-kV  line,  with  a  A-grounded-Y  transformer  at  the  load  end.  The  dc  injection  in 
this  case  is  in  the  load  circuit.  For  this  configuration,  the  initial  in-rush  transient  with  no  dc 
injection  is  large  and  its  time  constant  is  several  tens  of  seconds.  Figme  25  presents  the 
neutral  ciurent  measured  in  the  load  section  of  the  line. 
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Figure  25.  Neutral  Current  in  Transformer  Secondary  (Load)  for  Configuration  IB  with  No 
Dc  Injection. 

For  this  configuration,  an  attempt  was  made  to  determine  an  effective  relaxation  time 
constant.  It  was  observed  that  even  for  the  smallest  dc  current  injection  used  (5.0  A),  the 
system  relaxation  time  was  very  fast,  well  under  the  0.5-second  time  noted  for  the  previous 
case.  Evidently,  the  lowest  current  injection  had  already  driven  the  transformer  into 
saturation. 
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The  harmonic  content  of  the  12.47-kV  phase  1  current  and  voltage  is  shown  in 
Figures  26a  and  26b,  respectively,  as  a  function  of  the  dc  injection  level.  The  fact  that  the 
harmonic  content  of  the  current  jumps  by  about  1  order  of  magnitude  as  the  dc  current  level 
increases  from  0.0  to  5.0  A  indicates  that  the  core  saturation  occurs  somewhere  in  this  range. 
As  in  the  previous  case,  the  behavior  of  the  phase  1  neutral  voltage  is  not  affected  as  much  by 
the  dc  injection. 

Figure  27  presents  the  increase  in  reactive  power  on  phase  1  for  configuration  2. 
Evidently,  the  dc  injection  in  the  load  end  of  the  system  does  not  have  as  much  of  an  effect 
on  changing  the  reactive  power  demand  as  does  the  dc  injection  on  the  high-voltage  side  of 
the  transformer,  as  was  shown  in  Figure  23. 
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Figure  26.  Harmonic  Content  vs.  Dc  Cuirent  Level  for  Phase  1  of  Configuration  2. 
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8.0  CONCLUSIONS 


The  basic  measurements  made  in  this  experiment  can  be  used  to  answer  most  of  the 
questions  posed  in  Section  1.0  as  experimental  objectives.  Other  measurements  were  also 
made,  providing  additional  information.  Detailed  observations  of  the  test  results  are 
presented  in  Section  7.0. 

In  general,  if  the  dc  injection  is  limited  to  the  primary  side  of  the  transformer  as  in 
configuration  1 A  (Figure  2,  grounded  Y-Y)  the  majority  of  the  effects  are  seen  in  the  primary 
circuit.  However,  the  harmonic  distortion  at  lower-order  harmonics  (Figures  9b,  13,  14b,  and 
22)  will  be  passed  to  the  secondary  with  potentially  adverse  effects  on  the  connecting  load. 
Also,  with  sufficient  dc  injection,  the  generator  supplying  the  transformers  is  subjected  to  a 
very  high  reactive  power  demand  (Figure  16). 

For  cases  in  which  the  dc  injection  is  on  the  secondary  side  as  in  configurations  IB 
and  2  (Figure  3,  grounded  Y-Y,  and  Figure  4,  A-grounded-Y),  the  harmonic  content  of  the 
secondary  side  is  quite  high,  and  saturation  occurs  with  a  relatively  low  level  of  dc  injection 
(Figure  14).  However,  the  reactive  power  demand  on  the  primary  side  is  significantly  lower 
for  these  configurations  (Figure  27). 

Transformer  damage  is  dependent  on  the  dc  injection  level  and  direction  and  on  the 
core's  magnetic  history  and  thermal  performance.  During  this  experiment,  there  was  no 
apparent  damage  to  the  test  transformers  from  the  series  of  multiple,  10-15  and  26  second 
tests  at  various  dc  injection  levels.  The  thermal  mass  of  coolant  in  the  test  transformer  was  so 
large  that  several  minutes  of  injection  at  saturation  level  would  have  been  required  to  cause 
physical  damage;  hence,  there  would  be  no  damage  expected  during  a  10-15  second  test.  It 
should  be  remembered,  however,  that  MHD-EMP  is  expected  to  have  a  duration  of  several 
hundreds  of  seconds,  and  that  damage  to  transformers  due  to  such  heating  cannot  be  ruled 
out. 


While  the  transformer  response  is  certainly  nonlinear,  the  response  time  of  the 
transformer  can  be  evaluated  as  if  the  system  were  linear.  The  time  constant  determined  in 
this  way  is  found  to  be  much  lower  than  that  seen  in  transformer  tests  at  lower  dc  injection 
levels  on  large  power  transformers.  This  suggests  that  response  time  is  strongly  dependent  on 
injected  current  levels  and  that  the  effects  of  higher  levels  of  dc  injection  for  a  shorter 
duration  could  produce  very  high  reactive  power  demands  and  harmonic  distortions  in  a  few 
tenths  of  a  second. 
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Finally,  the  possibility  that  the  added  reactive  power  load  would  result  in  the  blowing 
of  fuses  on  the  primary  circuit  appears  to  be  a  concern  for  configuration  lA.  For  example, 
Figure  19  clearly  indicates  that  a  current  imbalance  results  form  the  injection  of  dc  in  the 
neutral  of  the  primary.  Assuming  a  normal  full-load  current  of  3.5  A  per  phase,  phase  1  with 
5.5-A  dc  injection  is  found  to  carry  about  5  A  for  the  duration  of  the  injection.  Hence,  the 
possibility  that  on  an  actual  power  system  the  primary  fuses  will  blow  with  sufficient  dc 
excitation  is  very  real.  In  addition  the  neutral  currents  are  quite  large  under  even  low-level  dc 
injection  (Figure  20)  and  could  result  in  a  failed  neutral  conductor  or  high  level  harmonic 
injection  into  equipment  via  the  neutral,  which  could  cause  equipment  failure.  This  is 
referred  to  as  the  "smoking  neutral"  phenomenon  in  situations  generating  high  harmonic 
content  in  the  neutral  conductor. 
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